Fluids of platelike particles near a hard wall 
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Fluids consisting of hard platelike particles near a hard wall are investigated using density functional 
theory. The density and orientational profiles as well as the surface tension and the excess coverage 
are determined and compared with those of a fiuid of rodlike particles. Even for low densities slight 
orientational packing effects are found for the platelet fluid due to larger intermolecular interactions 
between platelets as compared with those between rods. A net depletion of platelets near the wall is 
exhibited by the excess coverage, whereas a change of sign of the excess coverage of hard-rod fluids 
is found upon increasing the bulk density. 

PACS numbers: 61.20.-p, 61.30.Gd, 82.70.Dd 



I. INTRODUCTION 

While many theoretical studies have been devoted to 
the understanding of the behavior of elongated hard col- 
loidal particles near a hard wall (see, e.g., Refs. IT]- M), 
suspensions of disc shaped hard particles near a hard wall 
have not yet been investigated. One reason is that exper- 
imentally a corresponding hard platelet model system, 
i.e., consisting of particles with a short range repulsive 
interaction has been lacking until recently. From a the- 
oretical point of view the platelet fluid problem appears 
to be difficult because the Onsager approach of truncat- 
ing the corresponding equation of state at second order, 
which is valid for thin rodlike particles, cannot be jus- 
tified for platelets ^. Recently, preparation methods 
have been developed for new types of platelet suspen- 
sions, which may indeed serve as model systems of hard 
colloidal platelets |lG-fl2|. In the present article we use 
density functional theory (Sec. II) to study the positional 
and orientational order as well as the surface tension and 
excess coverage of thin hard platelets near a hard wall 
(Sec. III). Particularly, the density functional used here 
includes a third order density term which is not present 
in the Onsager second virial approximation. 



II. MODEL AND DENSITY FUNCTIONAL 
THEORY 

We consider an inhomogeneous fluid consisting of thin 
platelets of radius i? in a container of volume V. The 
platelets are taken to be hard discs without attractive 
interactions. The number density of the centers of mass 
of the platelets at a point r with an orientation oj — {9, cj)) 
of the normal of the platelets (see Fig. |l]) is denoted by 
/9(r, w). The equilibrium density profile of the inhomoge- 
neous liquid under the influence of an external potential 
V{r,uj) minimizes the grand potential functional |13| ]: 



FIG. 1. The polar angle 9 and the azimuthal angle of the 
normal u; of a thin platelet of radius R with its center of mass 
located at r = (0, 0, 0). 



n[p{r,Lo)] = / drdujp{r,uj) [fcijT (ln(47rA^p(r,cj)) - l) 

-^l + V{r,u;)]+Fe4p{r,Lo)], (2.1) 

where A is the thermal de Broglie wavelength and fi is 
the chemical potential. We express the excess free energy 
functional i^ex [p(r, i-u)] as an integral over all spatial and 
orientational degrees of freedom for the corresponding 
local free energy density fex{p = p{^,^)) of the fluid: 

Fex[p{r,uj)] = / dridLJidr2dLj2w{ri,r2,uji,uj2) 

xp{ri,uji)fea;{p{r2,uj2)), (2.2) 
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This functional incorporates fluid correlations via a 
weight function w{ri, r2,uji,W2) and bulk thermodynam- 



ics via a recently developed equation of state of thin 
platelets ^■. 



(2.4) 



where pf, — V^^ J dr J dto p{r,u!) is the number den- 
sity of the homogeneous and isotropic bulk fluid. 
We note that the resulting excess free energy func- 
tional i^es[p(r, w)] includes a cubic term of the den- 
sity, which is not present in the Onsager second virial 
approximation 0] 



fg°''Hp)^n'DL^pkBT 



(2.5) 



used in the description of thin rods of length L and diam- 
eter D near surfaces (see, e.g., Refs. [Q,|[ and Q). The 
necessity for including this higher order density term al- 
ready follows from comparing the ensuing bulk pressure 
(Eq. (2.4)) with corresponding simulation data MJ16| 
(see Fig. |]). Fluids consisting of thin platelets exhibit an 
isotropic and a nematic phase with no other liquid crys- 
talline phases, such as a columnar phase, observed [[l5[ . 
The isotropic-neniatic (IN) transition is first order with 
coexistence densities at pbiB? — 0.46 and pbNR^ = 0.5 

Minimization of il with respect to p{y, lo) leads to the 
following integral equation for the equilibrium density 
distribution: 

fcBTln[47rAV(r,w)] ^ p- V{v,uj) 
—Stt^R^ I dvidujip{ri,uJi)w{Y,Yi,uj,Uiii) 

^+^^'i?'[2p(r,c^)-fp(ri,c^i)] keT . (2.6) 

We have solved this equation numerically for a given 
chemical potential p and a given external field V{y^uj). 
The weight function is taken to be a function of the 
relative positions ri2 — ri — r2, and is normalized so 
that J dri2 (iw2 w;(ri2, Wi, W2) = 1- In this work the 
weight function is given by the Mayer function, except 
for the normalization. The Mayer function equals -1 if 
the platelets overlap and is zero otherwise. Two platelets, 
separated by a distance ri2, intersect if the inequality 



ri2 • [ijJi X oj2)\ < jR^sin-ff^ - (i"i2 •'^2)^ 



i?2sin722-(ri2-a;i)2, (2.7) 



with Ljj = {sm9jCOS(j)j,sin9jSin(l)j,cos9j), j — 1,2, 
is fulfilled. 712 is the angle between the normals Wi and 



has been derived and 
In practice, first the 



UJ2 of two platelets. Equation 

is amply documented in Rcf. 

weight function is calculated by testing Eq. (2.7) and is 

stored for all required values of (ri2,cji,cj2)- 







FIG. 2. Equation of state of thin platelets as obtained 
from Eq. (p.4| ) (solid line), and from computer simulations 
(squares |l5|] and diamonds u9). The dashed line follows 
from Eq. ^4[) by omitting the cubic term 



P?. 



Thereafter the integral equation (^) is solved using a 
Picard scheme with retardation. 



III. PLATELETS AND RODS NEAR A HARD 
WALL 

For model systems of hard particles near a hard wall 
at 2; = (see Fig. ^), apart from a possible surface freez- 
ing at high densities, nonuniformities of the density occur 
only in the z direction, so that p{r,uj) = p{z,9,(l)). Fig- 
ure displays the calculated platelet density profile for 
the bulk density ptR^ — 0.2. The calculations render 
p(z, 0, (j)) to be independent of the azimuthal angle for 
this density, i.e., there is no biaxial order emerging at 
the wall like for hard rods at high densities P-pf. For 
z < R orientations with large 9 values are forbidden so 
that the density profile has a discontinuity along the line 
Zmin = Rsiii9. There is a pronounced increase of the 
density near the surface. Moreover, slight packing effects 
are visible through oscillations at z « 2R. For compari- 
son Fig. ^ displays the density prof ile of t hin r ods near a 
hard wall calculated from Eqs. (2.1), ( ^.2D and (2.5). The 
bulk density has been fixed such that the second virial 
coef ficien t of the equation of state of thin platelets (see 
Eq. (2.4)) and of the equation of state of thin rods. 



Pb[l + -DL^pAkBT 



(rod) 
Pb 



are equal in units of R^ and ZJL^, 
Due to the presence of the wall p{z, 9, 

Z < Zmin = {L/2)C0S9. 



(3.1) 

respectively. 
) vanishes if 



y/R 2 




z/R 




x/R 




FIG. 3. The two systems under consideration consist of 
fluids of monodisperse thin platelets of radius R and thin 
rods of length L = 2R in contact with a planar hard wall at 
z — 0. Particles very close to the wall must adopt nearly a 
fully parallel alignment. Here the three principal directions 
are shown. 

Hence, in terms of the normals of the platelets and the 
normals of the rods along their main axis of symmetry 
the preferred orientation of rods near a hard wall is per- 
pendicular to the preferential orientation of the platelets. 
This means that in both cases the main body of the par- 
ticles tends to lie parallel to the wall (see Fig. 0). In 
agreement with earlier calculations [gj , for rods no pack- 
ing effects are visible because of the relatively smaller in- 
termolecular interactions between rods as compared with 
those between platelets. (For thin rods the intersection 
volume is pointlike whereas for discs it is like a line seg- 
ment.) 

We have tested successfully the accuracy of the numer- 
ical calculations by comparing our results for the profiles 
with the pressure sum rule 



Pb = ksT 



FIG. 4. Center of mass density profile p{z,6,cj))R of thin 
platelets near a hard wall for various angles 9 of the platelet 
normals relative to the surface normal. The density profile is 
independent of the azimuthal angle (j). At small distances z 
from the wall, large values of 9 are forbidden due to overlap. 
Therefore the profile is exactly zero if 2; < Zmin = -Rsin^. 



The bulk density is fixed to phR 



0.2. 



dO sinO p{z.rmn,d,4') ■ 



(3.2) 




and we found good agreement. Equation (3.2) is an ex- 
tension of the pressure sum rule of a hard-sphere fluid 
near a wall (see, e.g., Ref. iQ). The weighted density 
functional theory guarantees that the wall sum rule is 
satisfied, provided that the pressure enters through the 
corresponding bulk equation of state. 

A set of position-dependent order parameters quanti- 
fies the deviation from isotropy of the number density. 
For a uniaxial density profile the most general form of 
p{z, 6, (j)) is independent of (j) so that one can write: 



FIG. 5. Center of mass density profile p(^z^9^<\>)DI?' I\f^ 
of thin rods {L/D -^ 00) near a hard wall for the bulk 
density pbDL^ /\/32 — 0.2. The profile is exactly zero if 
z < Zmin ~ {L/2) COS 9 duc to Overlap. 



°° ■?/ 4- 1 



(3.3) 



;=o 



where P/(cos(6')) are the Legendre polynomials. The nor- 
malized, orientationally averaged number density profiles 
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FIG. 6. Normalized, orientationally averaged density pro- 
file niz) for platelets (lower solid curve) and rods (upper solid 
curve) in contact with a hard wall for the same bulk density as 
in Figs. W and H. The length L of the rods and the diameter of 
the platelets 2_R are taken to be the same. The dashed lines 
represent the results for an ideal gas of platelets and rods, 
respectively. The upper curves are shifted by 1. 

n{z) — 1-kQq{z)I p\, are displayed in Fig. || together 
with the results for non-interacting platelets and rods. 
Upon increasing z from the wall the averaged number 
densities increase and exhibit cusps at z = _R for platelets, 
and at 2 = L/2 for rods, respectively. For rods the 
maximum at the cusp is only about 10% above the bulk 
value which is essentially reached already for z = L. For 
platelets the maximum is more pronounced and slight 
packing effects are visible at larger values of z due to the 
relatively larger steric interactions between platelets. 

The position-dependent uniaxial, relative nematic or- 
der parameter s{z) = Q2{z)pb/{'^''^Qo{z)) is displayed 
in Fig. 0. At small values of z the value of the ne- 
matic order parameter reflects the geometric constraints. 
A platelet [rod] lying very closely to the wall must 
adopt nearly a fully parallel alignment, so that the or- 
der parameter reaches the limiting value 1 = P2(cos(0)), 
[—0.5 = P2(cos(7r/2))] there, whereas isotropic orienta- 
tion at large distances from the wall is characterized by 
s{z) = 0. Interestingly, the nematic order parameter for 
platelets has a minimum at z « 2i? due to a depletion of 
platelets parallel to the wall. In other words, a platelet, 
next to the platelet at 2 = i? with the rim touching the 
wall (see Fig. ||), is oriented rather parallel than perpen- 
dicular to the former. 

It is useful to consider not only local but also global 
properties of a liquid near a surface. There are two global 
quantities, which are of experimental interest and of in- 
terest for simulations: the excess coverage, which is 
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FIG. 7. Relative uniaxial nematic order parameter s{z) for 
platelets (solid curve) and rods (dashed curve) in contact with 
a hard wall for the same bulk density as in Figs. W and H. Posi- 
tive (negative) values of the nematic order parameter indicate 
that the platelets (rods) are preferentially aligned parallel to 
the wall. 

accessible by, e.g., gravimetric measurements; and the 
surface tension which, for example, is important for con- 
tact angles. The surface tension 7 is defined via 



fi[p(2,0,0)] = Fwb-|-75 



(3.4) 



where S is the surface area and lo\, — —pt (see Eq. (E^)) 
is the bulk grand-canonical potential density. The surface 
tension depends on the definition of what is denoted as 
the volume V [|l^. We have defined V as the volume 
of the container with its surface given by the position of 
the rim of the particles at closest approach so that the 
region < 2 < Zmin is part of V. Figure || displays the 
calculated surface tension together with the results for 
non-interacting platelets and rods. The steric interaction 
between the particles, which is more pronounced for the 
platelets, increases the surface tension with increasing 
density. The results for the rods are in agreement with 
those obtained in Refs. g| and |5|. 
The excess coverage 



r = 



00 

Pb / dz {n{z) - 1) 



(3.5) 



provides an important overall characterization of the den- 
sity profiles. Figure ^ summarizes the results for platelets 
and rods. Repulsive interactions between the platelets 
and the hard wall lead to a net depletion near the sur- 
face (r < 0). The excess coverage of rods exhibits a 
change of sign with increasing density in agreement with 
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FIG. 8. Surface tension of a fluid consisting of thin hard 
platelets of radius K (lower solid curve) and hard rods of 
length L and diameter D (upper solid curve, LjD -^ oo) 
near a hard wall. The dashed lines represent the results for 
an ideal gas of platelets and rods, respectively. The upper 
curves are shifted by 0.5 (dotted line). 



IV. SUMMARY 

We have applied a density functional theory to fluids 
consisting of thin hard platelets and rods in contact with 
a hard wall. Particles lying very close to the wall must 
adopt nearly a fully parallel alignment due to interac- 
tions with the wall (Fig. |3|). The probability of finding 
a particle touching the wall is increased compared with 
the bulk (Figs. ^ and ||). A comparison between the rod 
fluid and the platelet fluid exhibits slight orientational 
packing effects for the platelet fluid (Fig. |^) and that the 
increase of the surface tension with increasing density is 
more pronounced for platelets than for rods (Fig. |8|) due 
to larger intcrmolecular interactions between platelets as 
compared with those between rods. The calculated ex- 
cess coverage (Fig. ^ reveals a depletion of platelets close 
to the wall (Fig. O), whereas a change of sign of the excess 
coverage of hard-rod fluids indicate the onset of wetting. 
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FIG. 9. Excess coverage (see Eq. (B.5)) of a fluid consisting 
of thin platelets of radius R (lower solid curve) and rods of 
length L and diameter D (upper solid curve, LjD — * oo) 
near a hard wall. The dashed lines denote the corresponding 
results for an ideal gas of platelets and rods, respectively. The 
upper curves are shifted by 0.5. 
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